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Целью работы является построение математической модели, которая позволит проанализировать процесс 
загрязнения атмосферы. В качестве исходных данных взяты материалы государственных докладов, опубли- 
кованных в журнале «Экологический вестник Дона». Объектом исследования являются данные по выбросам 
загрязняющих веществ в атмосферу от автотранспорта, зафиксированные в Ростове-на-Дону и Ростовской 
области за период с 2001 по 2010 год. Математическое моделирование процесса загрязнения атмосферы 
осуществляется с помощью двух подходов: 1) описание с помощью математических моделей типа Гаузе; 
2) уточнение и прогнозирование процесса загрязнения атмосферы с помощью метода Прони. Используются 
программные пакеты МаЁСаа, Ма На. Результаты математического моделирования позволяют проводить 
анализ процесса загрязнения и выявлять его закономерности. 
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Введение. За последнее десятилетие наблюдается положительная тенденция в развитии эконо- 
мики. Выросли объёмы производства в машиностроительной сфере, в лёгкой и пищевой промыш- 
ленности, в строительстве и т. д. Как следствие, растёт благосостояние населения. Например, 
автопарк легковых автомобилей Ростовской области вырос за последнее десятилетие практически 
в два раза, что привело к росту выбросов в атмосферу вредных веществ. В настоящее время для 
уменьшения количества выбросов в окружающую среду применяются различные фильтры и 
устройства очистки. Тем не менее, объём выбросов вредных веществ в атмосферу продолжает 
расти, что следует из данных журнала «Экологический вестник Дона». 

Целью работы является построение математической модели, с помощью которой можно 
будет проанализировать процесс загрязнения воздуха, а также изучить закономерности этого 
процесса. 

Постановка задачи. В данной работе рассматриваются проблемы моделирования процесса за- 
грязнения атмосферы с использованием эколого-биологических моделей логистического типа и 
типа Гаузе, а также уточнение и прогнозирование развития процесса с помощью метода Прони. 
Для решения первой проблемы построены две модели: логистическая, которая не учитывает про- 
цесс очистки, и модель Лотки-Вольтерра, в которой процесс очистки уже учтён. 

Моделирование роста объёмов выбросов автотранспорта с помощью логистического 
уравнения. Атмосферный воздух является самой важной жизнеобеспечивающей природной сре- 
дой и представляет собой смесь газов и аэрозолей приземного слоя атмосферы, сложившуюся в 
ходе эволюции Земли. 

Используя многолетний мониторинг состояния биосферы, в том числе и атмосферы, про- 
водимый в журнале «Экологический вестник Дона» [1—6], можно сделать вывод, что наиболее 
существенным загрязнителем атмосферы в настоящее время является автотранспорт. 

Рост объёмов выбросов от автотранспорта наиболее точно описывается моделью роста 
популяции (логистическим уравнением) [7—10]. Логистическое уравнение запишем в виде: 


— =г(м 2) И (1) 
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где // — объём выбросов, Р — положительный параметр, г — скорость роста (увеличение объёмов 
выбросов), К — поддерживающая ёмкость среды (то есть, максимально возможный объём 
(концентрация) выбросов). 

Интегрируя уравнение (1), описывающее рост выбросов в зависимости от времени, получаем 

М=К № +Р. (2) 
М +(К-№)е” 

По статистическим данным, взятым из журналов «Экологический вестник Дона» [1—6], с помо- 

щью численных методов (численное дифференцирование, метод наименьших квадратов и т. п.) 

определяются параметры, входящие в уравнение (1). 

Метод наименьших квадратов (МНК) — метод оценки параметров модели на основании 
экспериментальных данных, содержащих случайные ошибки. В основе метода лежат следующие 
рассуждения: при замене точного (неизвестного) параметра модели приблизительным значением 
необходимо минимизировать разницу между экспериментальными данными и теоретическими 
(вычисленными при помощи предложенной модели). Это позволяет рассчитать параметры модели 
с помощью МНК с минимальной погрешностью. 

Мерой разницы в методе наименьших квадратов служит сумма квадратов отклонений дей- 
ствительных (экспериментальных) значений от теоретических. Выбираются такие значения пара- 
метров модели, при которых сумма квадратов разностей будет наименьшей: 

У =у.} > тт, 


7 
где У— теоретическое значение измеряемой величины, у— экспериментальное. 

При этом полученные с помощью метода наименьших квадратов параметры модели явля- 
ются наиболее вероятными [7]. 

Находим параметры К, ; №, В значения которых получаются равными К’`= 387,5, г = 0,6, 
№ = 8,1, Р= 186,3. 
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Рис. 1. График статистических данных (ломаная кривая, выделенная красным цветом) и решение логистического 
уравнения (выделено синим цветом) 


Логистическое уравнение достаточно точно описывает рост объёмов выбросов от авто- 
транспорта, что можно видеть из рис. 1. Степень отклонения полученной модели от фактических 
данных равна 2,2 %. 
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Математическое моделирование процесса загрязнения воздуха на основе модели типа 
Гаузе. Опишем модель взаимодействия загрязнения с окружающей средой на основе модели 
Лотки-Вольтерра [8—10]. 
Пусть ХС — объём (концентрация) загрязнений в момент времени & У — объём уловленных 
загрязнений в момент времени & 
ь =аХ - ВХУ (3) 
У =-аУ +сху 
где а, В, с, а — положительные постоянные, характеризующие скорость изменения выброшенных 
и уловленных загрязнений и взаимодействия между ними. 
Определяем коэффициенты а, В, с 4 методом наименьших квадратов (который описан 
выше): 
1) При помощи численного дифференцирования получаем переопределённую систему 
14,0 =-161с + 33488а, 
21,0 =208а - 33488, 
14,0 =229а - 40075, 
17,0 =-175с + 400754, 
—6, 0 = 243За - 466566, 
11,0 =-192с + 46656а, 
—5 =-203с + 48111а, 
22,0 =237а - 48111Ь, 
12 =-198с + 51282а, 
12,0 =259а -— 51282, 
32 =-210с +56910а, 
29,4 =271а 569106, 
7 =-242с + 72696, 8а, 
—35,4 = 300, 4а - 72696, 8Ь, 
—29 = -249с + 659854, 
—10,0 =265а - 659856, 
—2 =-220с + 56100а, 
6,0 =255а - 56100. 
2) Решаем переопределённую систему методом наименьших квадратов [7]. Получаем нормальную 
систему, которая имеет следующий вид: 
—1,2.10*а+2,6.10°6 =-1,4.10°, 
5,8 .10°а-1,2.108Ь = 10651,2, 
386888с —10.10’а = -9455, 
—10.107с+2,6.10°а =2,1.10°. 
3) Находим коэффициенты а, В, с а 


а:=0,4; 
Ь:=0,002; 
с:=-0,4; 

а =-0,001. 


Проверяем адекватность модели (рис. 2). 
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Рис. 2. Модель Лотки-Вольтерра 


Би -х+,-и 
>. хр +у? 


ворит о том, что модель вполне адекватна, особенно на начальном отрезке времени, но жела- 
тельно уточнение на последних шагах. 

Моделирование с использованием метода Прони. Для более точного исследования процес- 
са загрязнения применим метод Прони [11] — метод анализа коротких отрезков сигнала, осно- 
ванный на аппроксимации сигнала конечной суммой комплексных экспонент. Сигнал х[л] заменя- 
ется некоторой детерминированной экспоненциальной моделью: 


[п] = ХА ехр| (а, + 721, )(п-1)Т + 19, |, 


где Аки ах — амплитуда и коэффициент затухания (в С“. Ки 9+ — частота (в Гц) и начальная 
фаза (в рад) А-Й синусоиды. 








Получаем невязку по формуле ‚ Которая равна б =9%. Это го- 


Иначе Х[п]= Зульр ‚где №, =А, ехр(7-6,),2, = ехр (а, + 1211, )7Т ]. 


Задача заключается в поиске параметров экспоненциальной модели, аппроксимирующей 
р экспонентами 2р отсчётов действительных данных. 

Расчёты проводились с помощью пакета Ма аБ. Результаты вычислений представлены на 
рис. 3. Сумма экспонент достаточно точно отражает экспериментальные данные. Метод Прони 
имеет преимущества при исследовании апериодических сигналов по сравнению с другими мето- 


дами, а также может применяться для прогнозирования (на графике прогноз указан точками). 


Метод Прони 
т т т т т 





х м оби сибири, 














Рис. 3. Динамика выбросов загрязняющих веществ в атмосферу от автотранспорта в г. Ростов-на-Дону с 2001 по 2010 годы 


(на графике обозначены точками) и применение метода Прони (на графике ломаная кривая) 
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Заключение. В результате изучения процесса загрязнения атмосферы был проведён анализ 
данных. Построена модель роста объёмов выбросов от автотранспорта, которая наиболее точно 
описывается моделью роста популяции, но не учитывает процесс очистки. Также построена адек- 
ватная модель процесса загрязнения атмосферы типа Гаузе, в которой процесс очистки уже 
учтён. Произведено уточнение и прогнозирование процесса загрязнения атмосферы с помощью 
метода Прони. 
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Бот 2001 Ю 2010. ТРе таёпетайса/ ятиавоп о! {пе ай сотфат/пайоп ргосез$ [5 ретогтеа Егоидн о ар- 
ргоасйез: 1) аезспрйоп Бу тайетайса! тодеб оЁ баиге {уре; 2) птргоуетепЕ апа Югеса$Ипд оЁ {пе ай сощат- 
павоп ргосез$ Бу Ргопу5 те!оа; МаСаа, МаМаБ ргодгат раскадез аге изеа. Тве та{петайса/ хутианоп 
гезий а/ои/ {о апа/уге те сотщат!пайоп ргосез$ апа ю 'Аеп ву #5 геди/а Нез. 

Кеуигогаб: а!г ро/ийоп, те{петайса! тодбЕ!|, (09/58 с едиавоп, Ргопу'5 тетод. 
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